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We present the algebra for density operators projected to a topological band of a three-dimensional 
(3D) system. This algebra generalizes to 3D space the Girvin-MacDonald-Platzman algebra for the 
densities projected to the lowest Landau level in the case of the 2D fractional quantum Hall effect. 
We provide an example of a model on the cubic lattice in which the chiral symmetry guarantees a 
macroscopic number of zero-energy modes that form a perfectly flat band, and explicitly construct 
the algebra for the density operators projected onto this topological dispersionless band. The al- 
gebra of the projected density operators is related to the emergence of noncommutativity of the 
spatial coordinates of particles propagating in 3D, similarly to the noncommutativity of coordinates 
projected to the lowest Landau level in 2D. The noncommutativity in 3D is tied to a nonvanish- 
ing 9-texm associated to the integral over the 3D Brillouin zone of a Chern-Simons invariant in 
momentum-space. Finally, we find conditions on the density-density structure factors that lead to 
a gapped 3D fractional chiral topological insulator within Feynman's single-mode approximation. 



I. INTRODUCTION 

The integer quantum Hall effect (IQHE^ is the first 
known example of a fermionic phase of matter charac- 
terized by a topological index that is directly connected 
to a physical observable. The index in this case is the 
sum of the (first) Chern numbers obtained for each of 
the fully filled Landau bands, and the associated physi- 
cal observable is the Hall conduct ance.l^ll^ The fractional 
quantum Hall effect (FQHE) ^ results from the effects of 
electron-electron interactions when the Landau levels are 
partially filled with electrons, for certain rational filling 
fractions.'^ More examples of topological states of mat- 
ter that are comprised of noninteracting fermions have 
been discovered recently,'^'^ and have been classified ac- 
cording to discrete symmetries they respect or not, and 
the dimensionality of space in which the particles propa- 
gate. ^''^^ Such classification is sometimes referred to as 
the "periodic table" of topological insulators.'^ Among 
these states are topological ones associated with the 
presence of time-reversal symmetry (TRS) in two- as well 
as three-dimensional (3D) systems. It is natural to then 
question what the "fractional" version of these phases 
should be, and how they could be described. In partic- 
ular, it is interesting to ask what are the possible frac- 
tional topological phases of interacting fermionic systems 
in three spatial dimensions. 

One approach to capture universal physics arising from 
topologic al int eracting electron systems in (2-|-l)|i2ttlll 
and (3-1-1)21128] dimensions of space and time is via the 
parton construction: a fractional phase of electrons is 
obtained by constructing integer filled bands of "par- 
tons", which are then "glued" together by very strong 
gauge-mediated interactions so as to assemble together 
the physical electron. This approach is a generaliza- 
tion of theories that capture the universal physics of the 
FQHE, and yields, for instance, wavefunctions describing 



states which have fractional magneto-electric effects '^ZMl 
in the case of the topological insulators. The parton 
construction is one way to obtain an effective topologi- 
cal quantum field theory (TQFT) to describe fractional 
topological insulators. 

However, TQFTs do not capture the dynamics of the 
systems beyond their topological properties. As empha- 
sized by Haldane,^^ TQFTs are incomplete theories of the 
FQHE, for while they characterize the quantum numbers 
of the elementary excitations (topological defects), such 
as their charges and statistics, they do not contain any 
information about their energies. The information about 
the fundamental length scale in the FQHE, the magnetic 
length, is lost in its TQFT treatment. Recently Haldane 
in Ref s . and I5D1 has proposed a geometric description of 
the FQHE based on the algebra of the density operators 
projected to the lowest Landau level that was originally 
introduced by Girvin, MacDonald and Platzman (GMP) 
in Ref. [HI 

When projected to the lowest Landau level, the den- 
sity operators do not commute. However, the algebra 
closes in that the commutation of two density operators 
is proportional to a third one. Using this algebra, GMP 
were able to employ an approach that parallels that of 
Feynman and Bijl in their study of excitations in '^He.'^ 
Their approach allows to place a variational estimate on 
the excitation gap, if the static structure factor is known. 

The main objective of this work is to identify the den- 
sity algebra obeyed by electron operators in 3D topo- 
logical insulators, its relation to topological invariants, 
and its relevance to possible interaction-driven topolog- 
ical fractional phases in fermionic 3D systems. Armed 
with this algebra, one can forge ahead in trying to con- 
struct a dynamic theory of 3D fractional topological in- 
sulators that could perhaps parallel the solid understand- 
ing of the FQHE in 2D. In particular, the approach can 
suggest which types of interactions can give rise to in- 



2 



compressible gapped phases. 



II. 



PROJECTED ALGEBRA 



There are important symmetry considerations that 
need to be carefully taken into account when searching 
for interacting topological insulators in 3D. The FQHE 
descends from the IQHE when Landau or Chern bands 
are partially filled. In turn, the 2D IQHE is a stable class 
of states characterized by a Z index (symmetry class A 
in the terminology of Ref . [TB|) , which has no symmetry 
left out to be broken. If the logic is that we are also 
to start from a noninteracting topological insulator in 
3D when constructing the interacting fractional counter- 
part, we need to look at systems which are topologically 
nontrivial in 3D space. One possibility is to start with 
topological insulators. This has been the choice in most 
works so far. Here, instead, we shall start from systems 
that have chiral symmetry, but that lack TRS (symmetry 
class AIII in the terminology of Ref. [T^ . The rational 
for this choice is twofold. First, from experience working 
on strongly interacting 2D topological insulators, we 
have observed that TRS is easily broken in favor of mag- 
netized states due to the Stoner instability, which is en- 
hanced in bands with nonzero topological invariant .'^^lEIl 
Second, because the 3D chiral systems are characterized 
by a Z- valued topological invariant, it might keep a closer 
parallel to the FQHE. Indeed, we shall show that the 
density algebra for this 3D model does depend on this 
Z- valued topological invariant. 

The approach of GMP is ideally suited to the situ- 
ation where density operators are projected into a dis- 
persionless band (for example the lowest Landau level in 
the case of the FQHE). Here, we shall give a concrete 
lattice model with chiral symmetry that contains an ex- 
actly flat topological band, on which we construct the 
projected density operators. The resulting algebra will 
depend on the nonzero integral over the 3D Brillouin zone 
of a Chern-Simons action in momentum-space. For this 
lattice model, the average Berry curvature over the entire 
Brillouin zone is zero. Hence, the type of 3D fractional 
topological insulator that we discuss is qualitatively dif- 
ferent from the FQHE, where the average Berry curva- 
ture over the Brillouin zone is nonzero. The nature of 
the fractional states we discuss are intrinsically 3D, and 
not layered 2D (i.e., weak topological insulators). 

This paper is organized as follows. We show in Sec. |TT] 
how the position and density operators for noninteract- 
ing fermions, if projected onto the occupied bands of their 
insulating ground state, can generate a noncommutative 
algebra. Although this is done explicitly in 2D and 3D 
space, the method applies to any dimension of space. 
A minimal microscopic 3D noninteracting lattice model 
that realizes the conditions necessary to establish the 
noncommutative algebra of Sec.|IT]is presented in Sec. |III| 
The role of interactions is then discussed in Sec. IIVI 



We begin by recalling some elementary facts about the 
quantum motion of a spinless electron confined to move in 
the plane spanned by the orthonormal unit vectors and 
62 perpendicular to an applied uniform magnetic field 
B = B e^, whereby 63 = A 63. Its quantum dynamics 
is governed by the single-particle (Landau) Hamiltonian 
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2to, 



-A{R) 



B = VAA{r), (2.1) 



where the momentum = {P^ , ) and position Rj = 

(i?]^,i?2) operators obey the canonical commutation re- 
lation 



ihS„ 



(2.2) 



with jjLjV = 1,2. Hence, neither do the components of 
the covariant derivative in position space 



n 



h 



-A{R) 



(2.3a) 



nor do the components of the conserved guiding center 



commute, for 



and 



X :=R 



ni,H2 



-e, A n 



£2 



(2.3b) 



(2.4a) 



(2.4b) 



— \J hc/{eB) the magnetic length. 



with (.^ 

An orthonormal basis of energy eigenstates of the Lan- 
dau Hamiltonian (2.1 1 is made of the kets 



1 



\/n!m! 



(at)"(St)™|o), 



(2.5a) 



with 

St : 



V2h 



Hi + iH. 



B 



X^ + iX.^ 



(2.5b) 

and where n = 0, 1, 2, • • • labels the Landau levels with 
energy e„ = hujc {n + l/2) and m = 0, 1, 2, • • • , [($/$g)- 
1] labels the orbital angular momentum. Here, $ = AB 
is the magnetic flux threading the area A of the system, 
$Q = hc/e is the flux quantum, and lJc ~ eB/m^c is the 
cyclotron frequency. 

Defining the projector on the n-th Landau level 



(2.6) 
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one finds that the guiding center defined in Eq. (2.3b) equivalently that {/i,/2}p (''") = 1- In this case, it fol- 



is the position operator projected on any single Landau 
level 



lows that 



+i4 



V2 \ib - i6t 



From the projected coordinate algebra one can obtain 
a (projected) density algebra, by defining the projected 
density 



V2 " 



iS-i? 



ia^ — ia 



(2.7) 



p(r) :=P„e(r)P„, 
where the unprojected density operator is 



■V RV 



g{r) 



since V^aVn — and V^a'^Pn — 0, while V^bPn = b 
and Pnb^Pn — b^ ■ Thus, the position operators pro- 
jected to any given Landau level satisfy the noncommu- 
tative single-particle algebra (2.4b). 



Sir-R 



(2.11a) 



(2.11b) 



One can also naturally construct the guiding center oper- 
ators ( 2.3b ) from the projected density operators through 



This algebra is at the heart of both the IQHE and 
the FQHE. For example, it is intimately related to the 
quantized Hall conductivity a^. The Kubo formula for 
the contribution of the n-th Landau level (n = 0, 1, 2, • • • ) 
to the Hall conductivity is 



X^^ j d^rr^V^B{r)V^- (2.12) 

In momentum space, the projected normal ordered 
density operators 



(2.13) 



(2., 



where A is the area of the Hall droplet. This can be 
rewritten using Eq. (2.3a) as 

=2 



Ah 



2 r _ _ _ 

Y,[{n,m\RiKR2\n,m) ^ (1 ^ 2) 



in the lowest Landau level n = satisfy the commutation 
relations ^ 

[p{qi),PiQ2)] = -2i sin {q-^ A 93) ' piQi + ^2) 

(2.14a) 



X! X! f^"-' '^\^i'^n'R2\n, m) (1 O 2) Qj.^ ^j^g ^jj^j^ q£ g^g^jj wavevectors and q^, 



le 

Ah 



44 (q-i Aq2) •e3P(gi +92)- (2.14b) 



le 
Ah 



E 



h ' 



(2.9) 



where we used that A — 2tt^^£%. The role of the 
noncommutative position-operator algebra is apparent in 
the penultimate line. 

To quadratic order in the algebra of the projected 
position operators (2.3b I is maintained if a coordinate 
transformation — >■ f^{r), p, = 1,2, that varies on 
length scales larger than is area preserving. Indeed, 
we can then expand 



This algebra, the GMP algebra,'^^^^ plays two crucial 
roles. First, within the SMA approximation,^^ it dictates 
under what conditions interactions open a spectral gap 
between the many-body interacting ground state and its 
excitations upon lowering the chemical potential within 
the first Landau level. Second, it also dictates the uni- 
versal properties of the low-energy and long-distance dy- 
namics at the edge in an open geometry ! ""^^ 

The goal of the work presented in the remainder of 
this section is to generalize the noncommutative single- 
particle algebra (2.4b) to noninteracting many-body 



A(x),/2(x) 



(2.10a) 



fermionic Hamiltonians in 3D space. Before carrying out 
this program, let us motivate what it is to come by first 
presenting what would constitute a natural extension of 
the algebra in the QHE to 3D problems. 

First, instead of the commutator, consider the case 
where the 3-bracket of the 3D projected position opera- 
tors equals a C-number 



where the classical Poisson bracket is defined as 

ni. ,'dfl 9/2 



i^^ 



(2.15a) 



{/i,/2}p {r) 



Or,, Or, 



(2.10b) 



The condition for this coordinate transformation to lo- 
cally preserve area is that its Jacobian equals unity, or 



Xi , X2 , ^3 

where, following Nambu,'^we have defined the 3-bracket 

= [A,A^] A3 + [A2, 13] A, + [I3, li] I2. 

(2.15b) 
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The characteristic length scale I of the 3D noninteracting 
many-body Hamiltonian, not to be confused with the 
magnetic length 1^ of the 2D Landau Hamiltonian, is 
the signature of a spectral gap separating the ground 
state from the excited states. Similarly to the 2D case, 
where area preserving coordinate transformations left the 
commutation relations unchanged, we would like volume 
preserving transformations not to change the 3-bracket. 
Under generic transformations — /^(r), /i ~ 1,2,3, 
that vary on length scales larger than 



A (X), /2 (X), /2 (X) = {/i , /2, /gjj, (X) + 0{l% 



where the classical Nambu bracket is defined as 



2.16a) 



dr^ dr^ dr^ 



(r). (2.16b) 



The condition for this coordinate transformation to lo- 
cally preserve volume is that its Jacobian equals unity, 
or equivalently that {fi, /j, f^}^ {r) = 1. In this case, it 

follows that [/i(X),/2(X),/3(X)] ^i£^ + Oif-). 

Second, we posit that the 3D counterpart of the pro- 



jected density algebra (2.141 is 



[p(9i),p(92):P(g3)] = ©(^^ (91 Aq2)-93) KQ1+Q2 + Q3) 

(2.17) 

for some function Q with the Taylor expansion Q(x) cx x 
for small arguments x <C 1. 

The algebra defined by Eqs. plsl and (pTTl), if it 



can be realized by a 3D fermionic noninteracting many- 
body Hamiltonian, might then deliver two results. First, 
within the SMA approximation, it might dictate under 
what conditions interactions open a spectral gap between 
the many-body interacting ground state and its exci- 
tations upon lowering the chemical potential below the 
single-particle gap. Second, it might also dictate the uni- 
versal properties of the low-energy and long-distance dy- 
namics at the boundary in an open geometry. 

The key idea to realize the algebra defined by 



Eqs. (2.151 and (2.17) is to replace the effect of the mag- 



netic field in the Landau Hamiltonian by that of the pro- 
jection of suitable operators on a suitable subspace of 
the fermionic Fock space. The construction of this suit- 
able subspace presumes the existence of fermionic Bloch 
bands as occurs in condensed matter physics and assumes 
that a subset of these bands are fully occupied, while the 
complementary set are empty and separated from the 
filled subset by an energy gap. 

Now, carrying out this program for some Bloch bands 
will not yield immaculately the 3-brackets (2.151 and 



(2.17). It will yield these relations approximately in the 
long wave-length limit. The situation here is similar to 
the case of t he qua ntized Hall effect in flat Chern bands 
of 2D models.ESlSSl As discussed by Parameswaran, Roy, 
and Sondh i in Ref . l49l (see also Refs.l50land l5T|) . the alge- 
bra (2.14b) follows if the fluctuations in the Berry curva- 



ture over the Brillouin zone are neglected, or equivalently 



if the local curvature is approximated by its average over 
the entire Brillouin zone. Without this approximation, 
however, the projected position operator algebra will not 
be as simple as in Eq. (2.4b) and may instead be repre- 
sented as 



Xi,X2 



(2.18) 



where • • • stands for operators that appear as a result of 
the inhomogeneities in the Berry curvature. The central 
question is how to distinguish universal from nonuniver- 
sal contributions to the right-hand side of Eq. (2.18). To 



answer this question, we propose to consider the ground 



state expectation value 

quantized Hall conductivity, as seen in Eq. (2.9) 
show in Sec. Ill A] that 



that encodes the 
We 



1 



Xi, X2 



27ri 



Ch(^). 



(2.19) 



Here, Ch*-^-* is the first Chern number of the topologi- 
cal band that sust ains t he IQHE in the lattice, and will 
be defined in Eq. (2.46), while iVp is the total particle 
number and p is the average particle density. This sug- 
gests that the universal physical properties are captured 
by the C-number contribution to the right-hand site of 
Eq. ( [2l8| . 

As with the commutator (2.18), the 3-bracket (2.15) 



will also acquire extra terms in 3D space 



Xi,X2,X^ -i 



(2.20) 



We are thus lead to consider its normal ordered expec- 
tation value instead, which, as we show in Sec. |II A[ is 
given by 



(2.21) 



Here, the symbol CS^^^ stands for the 3D Chern-Simons 
invariant, which is indeed a quantized topological topo- 
logical invariant in systems with chiral symmetry, and 

will be defined in Eq. ( |2.44b| . 

We start by deriving Eq. (2.21) in Sec. 



II A for any 



Hamiltonian that is endowed with translation invariance, 
a spectral gap, and describes the motion of noninteract- 
ing fermions in flat Euclidean space R^. We then special- 



ize in Sec. II B to the case of massive noninteracting Dirac 
Hamiltonians for which this algebra simplifies greatly in 
the long-wave length limit. Finally, Sec. |II C] is devoted 
to the algebra of density operators in 3D lattice models, 
and will substantiate Eq. (2.17). 



Algebra of the position operators 



In Sec. |II A[ we study the consequences of the alge- 
bra of projected position operators in a situation where 
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FIG. 1: (Color online) Assumed spectral gap in the single- 
particle energy spectrum. Here, /j.^ denotes the chemical po- 
tential and the insulating noninteracting many-body ground 
state |<E>) is obtained by filling all the states in the A'^ bands 
below the spectral gap. 



fermions fully occupy the lowest N bands of a total of N 
bands of a noninteracting translational invariant (Bloch) 
Haniiltonian (see Fig. [l]) . We assume that the N bands 
are separated by an energy gap from the higher bands 
and that the chemical potential lies in this spectral gap. 
We denote with V~ the projection operator on the states 
in these gapped lower bands. In analogy to the guiding 
center coordinates (2.12) in the IQHE, we can then define 
the operators 



X,. 



V r T'- 



g{r)r 



N 



(2.22) 



in d-dimensional space for any fi = 1, ■ ■ ■ , d, where g{r) 
measures the fermionic density at position r. The engi- 
neering dimension of the operators is that of length. 
Their expectation value for any single-particle state is 
nothing but the expectation value of the /it-component 
of the projected position operator in this single-particle 
state. Hence, we interpret this operator as the fi- 
component of the projected many-body position opera- 
tor. When multiplied by the electric charge, it is nothing 
but the projected polarization operator. 

In the following, we will compute the expectation val- 
ues of the commutator and the 3-bracket of the operators 
X^ — I, ■ ■ ■ , d) in the insulating noninteracting many- 
body ground state |$) obtained by filling all the states 
in the N lower bands. For that, we will use the second 
quantized formalism, within which the density operator 



g{r) := V^i(r)V^„(r) 



(2.23) 



is expressed in terms of the second quantized single- 
particle operator tpl^ (r) that creates a fermion at position 
r and in the orbital labeled by a, where summation over 
the a = 1, ■ • • , N orbital degrees of freedom is implied. 
We arc using the normalization convention 



{V^aM>l.(r-')}='5„,o"5(r-r') 



(2.24) 



for the anticommutator between creation operators with 
orbital a = 1, ■ ■ ■ , N and position r e M'' and annihila- 



tion operators with orbital a' = 1, • • • , N and position 
r' € R"^. 

We denote with ip^ik) the Fourier transforms of the 
second quantized operators V'a('")i « — 1, • • • ,N. The 
projection V~ is most conveniently carried out after 
rewriting ip^i^) '^'^^ ^ unitary transformation 

^^ik)^ui^Hk)Xaik) (2.25) 

in terms of the operator Xai^^) that creates a state at 
momentum k e M"^ in the band labeled by a, where sum- 
mation over the a ~ 1, • • • , N bands is implied. For any 
given momentum k G M.'^ and band a = 1, • • • ,N, the 
vector with the N components u'^\k) labeled by the 
orbitals a — 1, • • • ,N is nothing but the eigenstate of 
the Bloch Hamiltonian. These Bloch states obey the or- 
thonormality condition 



N 



a,b=l,---,N. (2.26) 



Then, the projection V~ amounts to restricting the sum- 
mation over the bands a = 1 , ■ • • , N that is implicit in 
Eq. (2.25 1 to the subset of the a = 1, • • • ,N lowest bands. 



Hence, the projected density operator takes the form 
V,iir,P, = jly, I i'kf, ^^^^^ 

xxi(l=)"?"(fc)"?''('=')la.(«:'). 

where the symbol " is used to indicate that the implied 
sums over the repeated band indices a, a' run only over 
the lower band indices a, a' = 1, • • • , A^ on which the 

operator is projected. 

We can now use the form (2.27) of the projected den- 



sity operator to rewrite the /i-component of the projected 
position operator ( 2.22 ) as a single integral in momentum 
space,!^ 



X,^ I d''kxlik)xf{k)xi{k), 



(2.28a) 



where = 1, 
dinates 



, d, we have introduced the N x N coor- 



xf(fe)^il?f(fe):=i 6'%+Al\k) 



(2.28b) 



and a summation is implied over the repeated occupied 
bands a, 6 = 1, • • • ,N {d^ is a short-hand notation for the 
partial derivative with respect to fc''). These coordinates 
are nothing but the N x N matrix elements of the non- 
Abelian covariant derivative 



xJk) = iDJk) 



(2.28c) 



for any k £ M."^ in momentum space, since the non- 
Abelian gauge field A^{k) with the matrix elements 



Af{k) 



u^:>{k) (fc). 



(2.28d) 
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where fi = 1, ■ • • ,d and d,b = 1, • • ■ , N, realizes a non- 
Abelian Berry connection. 

The x^{k) are noncommutative if the non-Abehan 
gauge field has the nonvanishing field strength 

F^,{k) [D^{k),D,{k)] 

= {d^A,~d,A^) {k)+[A^{k),A,{k)] 

(2.29) 

for some = 1, ■ ■ • , d and some k € in momentum 
space. 

We are now going to construct a family of expecta- 
tion values of products of the projected position opera- 
tor (2.28a) in the noninteracting many-body ground state 

For that purpose, it is useful to relate the noninter- 
acting many-body ground state expectation value of any 
projected many-body operator O that is bilinear in the 
occupied band operators to a single-particle expectation 
value of its first-quantized kernel O. In other words, we 
consider the second-quantized operator 

d:=Jd''k j d'k' xl{k)0^^\k,k')xAk') (2.30a) 

where summation is implied over the repeated occupied 
band indices and, for any pair of occupied bands a, a' = 
1, • • • , iV, the first-quantized operator 0"'^{k, k') is some 
differential operator with respect to the wave number 
k' eR'^. It then follows that 



O 



= / d^'fc / d'^k' 



X ( <i>, fc,a 



O^^ {k,k') 



(2.30b) 



On the right-hand side, we have introduced the single- 
particle ket |$,fc,a) defined by removing from the non- 
interacting many-body ground state |$) one fermion 
with wave number k S and occupied band index 
5, = 1, • • • whereby these single-particle states are 
normalized according to the convention 



(<i>,fc,a|<i>,fc',a') =<5a,5,(5(fe-fc') 



(2.30c) 



for any a, a' = 1, 
normalized as 



, iV and fc, fc' e M'' and |$) would be 



($1$) = 1 



(2.30d) 



if periodic boundary conditions were imposed in a box 
of linear size L before the thermodynamic limit L oo 
was taken. Furthermore, we shall assume that 



0^''{k,k') = S{k - k')0^''{k') 



(2.31a) 



decomposes additively into two pieces for any a,b — 
1, • • • , N and k,k' G M''. The first piece is diagonal in the 
band index but may contain differential operators with 



respect to the wave number k' G M.'^. The second piece 
need not be diagonal in the band index but acts multi- 
plicatively (like a C number) on the band operators. The 
covariant derivative D^ik') from Eq. (2.28c) when mul- 
tiplied from the left by 5{k — k') satishes this additive 
decomposition. If so, we conclude that 



O 



= y d'^fc ($,fc,a|e''^"(fc)| $,fc,s) 
=AAx Jd'^ktiOik) (^-^^^^ 



We have introduced the single-particle wavefunction nor- 
malization constant 



Af := fc,a|$, k,d) 



(2.31c) 



that is independent of the momentum A; e M*^ and of the 
occupied band a = 1, • • • , N. This constant is the inverse 
volume element in momentum space, i.e., it would scale 
like L'' as L — ^ oo if periodic boundary conditions were 
imposed in a box of linear size L before the thermody- 
namic limit i — > oo was taken. We have also introduced 
the trace tr over the occupied band index and the dimen- 
sionless trace 



Tr^ (•••)= X y d'^fc tr (• • • ) 



(2.31d) 



over the single-particle Hilbert space on which O acts. 

The first member of the family of expectation values is 
that of the projected position operator for any fi — 
1, • • • ,d given by 



X, 



iTr.i^^. 



(2.32) 



The second member in this family is constructed from 
the second-order polynomials 



X ,, , X. 

fj, ' I 



for any pair fj.jV = 1, ■ ■ • 
Eq. (2.28a), one verifies that 



(2.33) 

d. After making use of 



d''kxl{k)F^^t{k)xi{k) (2.34) 



is a normal-ordered bilinear form in the x operators for 
any pair ii,v — 1, ■ • • ,d and with the summation con- 
vention implied over occupied projected band indices 
d,b — I,-- - ,N. Consequently, 



$ [X^,X,\ ^) = -Tt,F^^, (2.35) 

for any pair fi,!^ = I,-- - ,d. Observe that the non- 
Abelian contribution to F^^ defined in Eq. (2.29) drops 
out from 



tri^^,(fe)=tr {d^,A,-d,A^) (fc) 



(2.36) 
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for any fc G M'' and any pair fj,,i/ = 1, ■ ■ ■ ,d. 

The third member in this family is constructed from 
the normal-ordered third-order polynomials 



Then, Eq. (2.401 is invariant under those local U{N) 
gauge transformations (2.391 of the single-particle 



(2.37) 



Hilbert space that leave the boundary conditions un- 
changed. 



In contrast, Eq. (2.35) can be rewritten 



for any given triplet /x, i^, p — I,-- - ,d, written as the 



quantum 3-bracket defined in Eq. ( 2.15b I. Normal order- 



ing is defined by moving all creation operators to the left 
of the annihilation operators as if these operators were 
mere anticommuting numbers and is indicated by sand- 
wiching the expression to be normal ordered by colons. 
After making use of Eq. (2.28a) and Eq. (2.34), one ver- 
ifies that 



1 



J d^ktrF^,{k) (2.42a) 



for any = 1, - ■ ■ ,d. Here, we have multiplied the ex- 
pectation value by the inverse of total number of particles 
iVp to obtain an intensive quantity. Then, the ratio of the 
number of particles N and the single-particle wavefunc- 



-e"^^Tr,(i^,)Tr,(F^^) 
-iTr, [D^,D,,Dp] 



tion normalization constant Af defined in Eq. (2.31c 
nothing but the average particle density 



(2.38) 

for any given triplet fi, v, p = 1, ■ ■ ■ ,d where a summation 
convention is implied for the repeated indices I, J^K = 
H, V, p and we have introduced the fully antisymmetric 
Levi-Civita symbol with e^'^P = +1. 

One verifies that both the left-hand sides and the right- 
hand sides of each of Eq. Eq. (Essl), and Eq. 



p := 



(2.42b) 



Equation (2.42a) is invariant under the local U{N) gauge 
transformations (2.39) of the single-particle Hilbert 



space, even if they change the boundary conditions. 

Finally, the last term on the right-hand side of 
Eq. (2.381 can be rewritten as, following Appendix [A] 



are invariant under the local U (N) gauge transformation 
defined by 

xHk)^xHk)G\k), x(fc)^G(fc)x (fe), 



Tr, [D^,D„D^] 
3 



JJK 



A{k) G{k)A{k)G\k) - {dM){k) g1'(A;), 



(2.39) 



d'^fe tr ( AjFjK - -A,AjA„ j (fc), 

(2.43) 



where G{k) is any unitary N x N matrix for any k e 
M'' and matrix multiplication is implied in Eq. (2.39) 



where L J, K ^ p,v,p. As with^Eq. ( |2.40| ), Eq. ( |2.43[ ) 
only holds, if the local U{N) gauge transforma- 



with the operator- valued column- vectors and row- vectors 
X (fc) and x^(fc) that have the components Xhi^) ^^'^ 
Xa(fc), a = 1, • • • , iV, respectively. 



tions (2.39) of the single-particle Hilbert space are re- 



and (2.38) on the one hand, and Eq. (2.351 on the other 



stricted to those that leave the boundary conditions un- 
changed. 

Now, in order to relate the right-hand sides of 



T here i s a qualitative difference b etwee n Eqs. ( |2.32[ ) Eqs. (|2.40[), (|2.42a|, and (|2.43|) to topological invariants, 



hand. The right-hand side of Eq. ( 2.35 1 solely depends on 
the single-particle operator F^^ with p^v — 1, ■ • • ,d, an 
operator that acts multiplicatively on the single-particle 
basis |3>,fc) as far as the momentum dependence is con- 
cerned. This is not so for Eqs. ( |2.32[ ) and \2M\ . Their 
right-hand sides involve single-particle traces over differ- 
ential operators with respect to the components with 
p = 1, - ■ ■ , d in the single-particle basis |$, fe). 

The expectation value (2.32) consists of two contribu- 
tions 



we consider periodic boundary conditions in momentum 
space, instead of the open boundary conditions that were 
imposed up to this point. It follows that the domain of in- 
tegration in momentum space is a d-dimensional torus T'^ 
with the volume (27r)'^ in momentum space, we identify in 
the right-hand sides of Eqs. (|2.40|), (|2.42|), and (|2.43|) the 



X,. 



j d'ktrA^ik), (2.40) 



the first of which is ill-defined and refiects the arbi- 
trary choice of origin of the coordinate system. The 
second term, however, corresponds to the physical 
polarization.EH Here, we will fix the origin such that the 
derivative becomes a traceless operator 



topological invariants associated with the single-particle 
states that span the iV-tuple of bands. We define the d 
ID Chern-Simons invariants in rf-dimensional momentum 
space as 

CSl'^ -.^i-^tvA^, p=l,---,d, (2.44a) 

and the d{d—l){d—2) /6 3D Chern-Simons invariants in d- 
dimensional momentum space for any choice oi p^v, p — 
1, • • • , d as 



iTr^ 



d 



JJK f jrfL. 



0. 



(2.41) 



(2.44b) 
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where I, J, K = fi, p. The integral on the right-hand 
side of Eq. ( |2.44a[ ) and Eq. ( |2.44bp is quantized to half- 
integer values if the single-particle Hamiltonian obeys the 
chiral symmetry and the domain of integration is that of 
a d-dimensional torus T'* with the volume (Stt)'^. The 
ID and 3D Chern-Simons invariants in d-dimensional mo- 
mentum space carry the engineering dimensions of length 
raised to the powers (1— rf) and (3— d), respectively. They 
are thus dimensionless if and only if c? — 1 and d ~ 3, 
respectively. 



The Chern-Simons invariants (2.44a) and (2.44b) are 



only well-defined modulo int eger v alues, since the local 
U{N) gauge transformations (2.39) of the single-particle 



Hilbert space can change them by their winding numbers, 
that are 

i/^^tr Gta^G, = (2.45a) 



and 



6 J (2.)^ t^^' ^^^^ ^'^^ ' 

(2.45b) 

with I,J,K = fjL, V, p and p,,v, p = 1, - ■ ■ , d, respectively. 

In contrast, the d{d— l)/2 first Chern numbers defined 
in d-dimensional momentum space as 



ChW .^i / ^ 



tr F„„, p.,v^l,--- ,d, 

(2.46) 



can only take integer values if the domain of integra- 
tion is that of a c?-dimensional torus T"^ with the volume 
(271)** in momentum space,!^ irrespective of whether or 
not the single-particle Hamiltonian obeys the chiral sym- 
metry. However, when c hiral s ymmetry holds, the ID 
Chern-Simons invariants (2.44a) are quantized.^^ There- 



fore, derivatives of these quantities vanish, which in turn 
implies that all first Chern numbers (2.461 vanish. Fur- 



thermore, the first Chern numbers (2.461 are invariant 



under the local U{N) gauge transformations (2.39). The 
first Chern numbers defined in d-dimensional momentum 
space carry the engineering dimensions of (2 — d). 

Let us explicitly state our main result for the fully 
antisy mmetr ic product of three projected position oper- 
ators (2.28a) in the noninteracting many-body ground 



state for the case of d = 3 dimensions 



1 



(27r)2i 



^ X e^^'-P CSfi) ChW + 3CS(3) 



(2.47) 



[The single-particle wavefunction normalization constant 
A/" was defined in Eq. ( 2.31c ). The average particle den- 



sity p was defined in Eq. (2.42b).] Again, it should be 



noted that the right-hand side of Eq. (2.47) is entirely 



determined by its quantized topological numbers if the 
single-particle Hamiltonian obeys the chiral symmetry, 
periodic boundary conditions hold, and if the equality is 
understood modulo contributions from local U{N) gauge 
transformations (2.39) with non- vanishing winding num- 



bers, in which case 



Xi,X2, 



$) = -^^CS(^) {2.4 



P 



owing to Ch^,^-* — for any v, p 

We close Sec. HA by noting that it is natural to gen- 
eralize our discussion to d-brackets of projected position 
operators in d dimensional space defined by 



1, 



X X,. 



(2.49) 



where Pi, - ■ ■ , — 1, ■ ' • and e^i "^<i is the fully an- 
tisymmetric pseudotensor (Levi-Civita symbol) of rank 
d. Since Eq. (2.49) is antisymmetric in any two in- 



dices, one can infer from Eq. (2.34) that the noninter- 



acting many-body ground-state expectation value of the 
d-bracket (2.49) is expressed solely in terms of the Berry 



. ^^(fc) whenever d is even. By Wick's the- 



curvature F,, 

orem, this expectation value involves a linear superpo- 
sition of many-body contributions ranging from a d/2 
many-body contribution all the way down to a one-body 
contribution. The one-body contribution is the d/2-th 
Chern number which is proportional to 



pMlM2---Md-iA'd 



d-^fctr {F^^^^ X 



X F, 



(2.50) 



In contrast, for any odd dimension of space d, the nonin- 
teracting many-body ground-state expectation value of 
the d-bracket for the projected position operators in- 
volves Chern-Simons forms as was the case for c? = 1, 3 
as we have already explored. 



B. 



Massive Dirac fermions and volume-preserving 
diffeomorphisms 



In Sec. |II A[ we have related the ground state expecta- 
tion values of the commutator and of the 3-bracket of pro- 
jected position operators to quantized topological num- 
bers, namely the Chern numbers and Chern-Simons in- 
variants. By contrast, we have recalled in Eq. (2.4b) that 



a Landau level has the special property that the commu- 
tator of projected position operators itself is nothing but 
an imaginary number 



= - F, 



(2.51) 



where — 0,1. In other words, the Berry curvature 
is constant in a Landau level. The algebra ( 2.51| ) entails 
the so-called algebra obewd by the projected density 
operators in a Landau level.'^ 
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Here, we are going to show that the same is true for 
massive Dirac electrons in 2D, if the Umit of small mo- 
menta fc — 7> is considered. We then extend the discus- 
sion to massive Dirac electrons in 3D, where we consider 
the 3-bracket of projected position operators in the same 
limit of small momenta. 

In 2D Euclidean flat space, a single flavor of Dirac 
fermions with mass m and in the fundamental represen- 
tation of the Lorentz group is governed by the single- 
particle Hamiltonian in momentum space 



'^2D(fe) + ^2(^2 + Tna^ 



(2.52) 



As usual, we use Uq for the 2x2 unit matrix, while a^^ 
(72, and (Tg are the three Pauli matrices. 

This Hamiltonian supports two bands with the nonde- 
generate energy eigenvalues 



and the Berry curvatures 



m 



2[e(±)(fc)]' 



(2.53) 



(2.54) 



for II, v — 1,2. Upon projection to the lowe r band 
£^~^( fc), t he position operators defined by Eqs. (2.28b) 
and (2.54) satisfy the noncommutative algebra 



,,]{k)=-Fj,-Hk) 



fj,: 



^ sgnm -I- O(fc^) 



(2.55a) 



for II, v = 1,2. The Dirac counterpart to the magnetic 
length in the QHE is here 



1 



(2.55b) 



As announced, the algebra ( |2.55 ) reproduces the alge- 
bra (2.51 ) in the limit fc — > 0. The first Chern number of 



the lower band is given by 



1 

2n 

sgnm 
2 



(2.56) 



We will now illustrate a geometric consequence of the 
algebra (2.55) of the projected position operators. For 



that purpose, consider a pair of smooth functions fi (r) 
and /2 (r) that vary on length scales larger than t^y. With 
the help of the Taylor expansions of fi (r) and /2 (r) 
around the origin, we can define two single-particle oper- 
ators fi (x) and /2 (a;) by replacing r in the Taylor expan- 
sions with the projected position operator x = (a;„). As 



a consequence of the algebra (2.55), these single-particle 



operators satisfy in the limit of long-wavelength 

[A , h] i^) = sgn m {A , A }p (x) + 0(^d), (2-57) 



where the classical Poisson bracket was defined in 
Eq. (2. 10a I. The commutator (2.571 is invariant under 



area-preserving diffeomorphisms of the Euclidean plane, 
just as the Poisson bracket is, to leading order in an ex- 
pansion in powers of fp. The algebra obeyed by the set 
of diffeomorphisms of the Euclidean plane that leave the 
Poisson bracket invariant realizes the so-called classical 
algebra. Thus, Eq. (2.57) draws the connection to 
the quantum version of the algebra, which is the 
algebra (see Refs. EHl |39l HQl and IH]) ob eyed by 
the projected density operators in a Landau level. '^i'^^Hlll 
A manifestation of this correspondence is found in the 
nondissipative Hall viscosity, which can be viewed as the 
response function of the quantum fluid to an infinites- 
imal area-preserving deformation.ISfil In turn, an incom- 
pressible 2D classical fluid may be described in terms of 
a one-form gauge field, as appears in the Chern-Simons 
theory relevant to the quantum HaU effect (QHE).l^^ltSl 
In 3D Euclidean flat space, a single flavor of Dirac 
fermions with mass m and in the fundamental represen- 
tation of the Lorentz group is governed by the single- 
particle Hamiltonian in momentum space 



^3D(fc) — ~ i"^/^"^; 



5, (2.58a) 
where we have defined the Hermitian 4x4 matrices 







/3 



To 

-ar 



75 







(2.58b) 

Observe that this Hamiltonian has the chiral symmetry 



75'^3D(fc)75 = -'HMk) 



(2.59) 



for all fc G M'^. The spectrum of Hamiltonian (2.58a 
consists of two doubly degenerate bands with the energy 
eigenvalues 



e(=^)(fc) ±v/fc2 +to2 (2.60) 
and the non-Abelian Berry field strengths 



Fl^\k) = ±Ule^,^^^+0{\k\) 



(2.61) 



for II, V = 1,2,3, where 7^ = (— U]^, (T2, Ug). Conse- 
quently, upon projection to the lower band e'-^''(fc), the 
position operators defined by Eqs. (2.28b I and (2.61 ) sat- 
isfy the noncommutative algebra 



[x^,xMk) = -dle^,^-t^ + 0{\k\) 



(2.62) 



for iJL,v = 1, 2, 3. The definition (2.55b) of has carried 
over. Since the non-Abelian Berry field strength (2.61) 



is traceless, all first Chern numbers (2.56) vanish as ex- 



pected for a system with chiral symmetry. In contrast, 
the 3-bracket is again proportional to a Dirac-like oper- 
ator 



[Xi,X2,X^(k) 



^D7%(fc)+0(|fe|) 

.3^3 



^d7^5^ 



0{\k\). 



(2.63) 
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As for the case of Dirac electrons in 2D, we can as- 



sociate a geometrical meaning to the algebra (2.631 of 



the projected position operators. To that end, consider 
three smooth functions /^(r), f2{'f), and f^{r) that vary 
on length scales larger than Again, with the help of 
the Taylor expansions of these functions, we define three 
single-particle operators by replacing r in the Taylor ex- 
pansions with the projected position operator x = {x^). 



As a consequence of the algebra (2.63), these operators 



satisfy in the limit of long wavelength 

tr[/i,/2,/3]N-i3\/2^^ {/i,/2,/3}n {x)+0{ti), 

(2.64) 

where tr denotes the trace over the matrix degrees of 
freedom and the classical Nambu bracket was defined 
in Eq. (2.16b). The Nambu bracket is invariant under 



volume-preserving diffeomorphisms of Euclidean space 
just as the Poisson bracket is invariant under area- 
preserving diffeomorphims in 20.*^ In the description of 
ideal 3D classical fluids a two-form gauge field naturally 
arises as a consequence of volume preserving diffeomor- 
phisms.'^ Such a two-form gauge field also appears in 
the 3D BF theory that is believed to be rele vant to 3D 
topological insulators-fSS Since the 3-bracket (2.64) is in- 



variant under volume-preserving diffeomorphisms of 3D 
Euclidean space, to leading order in an expansion in pow- 



ers of Eq. (2.64) draws a connection to a quantum 



algebra that generalizes the classical algebra obeyed by 
volume-preserving diffeomorphisms . 



C. Algebra of the density operators 

Until now, we have considered the algebra obeyed by 
the projected position operator assuming translation in- 
variance in Euclidean flat spaces. What are the conse- 
quences of the noncommutative coordinates in 3D space 
for the algebra of the projected density operators? 

To answer this question, we resort to a tight-binding 
model defined on a lattice A with a Brillouin zone BZ, 
and on which we impose periodic boundary conditions. 
We assume, without loss of generality, that the lattice is 
three dimensional. In this spirit, we turn our attention to 
the single-particle electronic density defined on a given 
site r of a lattice A as 



N 



Qr — \r,a){r, 



(2.65a) 



Q = l 



where a = 1, • • • ,N labels degrees of freedom on every 
lattice site, e.g., spin or orbitals. These operators obey 
the closed algebra 



dr.dr^ = Sr^.r^Qr, (2.65b) 

owing to the orthonormality of the single-particle states 
(ri,Qi|r2,a2) ^ S^^^^J^^^^^ (2.65c) 



for any pair of sites and from the lattice A and for 
any pair of orbitals ai,a2 = 1, • ■ ■ N. As a consequence, 
these operators commute pairwise. 

The Fourier transform of in terms of the orthonor- 
mal Bloch states \k,a) labeled by the wave number k 
from the BZ and orbital index a — 1, - ■ ■ , N reads 



N 



(2.66a) 

fceBZ Q=i 

for any q G BZ. These operators obey the closed algebra 

Qq,eq^=Qq,+q^^ (2.66b) 

owing to the orthonormality of the single-particle states 

(gi,ai|q2,"2> = ^q^,q^ ^a^a^ (2.66c) 

for any pair of wave number Qi and ^2 from the BZ and 
for any pair of orbitals ai, a2 — I, ■ ■ ■ N . As a conse- 
quence, these operators commute pairwise. 

Consider now a basis transformation in the a degrees 
of freedom for every k G BZ that is parametrized by 
the N X N complex- valued numbers with a,b = 

I,--- ,iV, i.e., 



N 



Y.utl\k,a), 6=1,..., TV. (2.67) 



The ket \uf.,h) labeled by fe € BZ for any given b = 
1, . . . ,N should be thought of as Bloch state of the b- 
th band of a single-particle Hamiltonian. This Hamilto- 
nian shares the translational symmetry of A and periodic 
boundary conditions are imposed. For any q G BZ, we 
define the density operator projected on a single (nonde- 
generate) band b by 



N 



EE 

feeBZ Q=i 



(b)* (b) 

"fc,Q '^k+q,a 



(2.68) 



The projected operators with q G BZ do not obey any- 
more the algebra (2.66b). In the long- wavelength limit 



9i ; 92 ^ 1 (the lattice spacing of A is set to unity) , their 
commutation relation is "^^^^^ 



Pq,:Pq, 



1ll2 E ^Mi^.fe 



feeBZ 



,{b) 



k+q^+q^ 



(2.69a) 



to leading order in an expansion in powers of the compo- 
nents of qi and ^21 where 



and 



N 



(2.69b) 



(2.69c) 
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for fj,,^ = 1, 2, 3 are the Abelian Berry curvature and 
the Abelian Berry connection, respectively, and is un- 
derstood as the derivative with respect to the momen- 



tum component k'^. The long- wavelength algebra (2.69a 



closes only if F^^ ^, is independent oik, in which case 
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— (qiAq2)-Chpg^+q^ (2.70a) 



to leading order in an expansion in powers of the compo- 
nents of Qi and where 



Ch^ 



27ri ef""^ 



V F 



(2.70b) 



fceBZ 



with A = 1,2,3 are the components of the vector Ch 
made of the three first Chern numbers characterizing any 
nondegenerate band in 3D space.^^ (A summation con- 
vention is implied for the repeated indices fi,i> = I, 2, 3.) 



In the thermodynamic limit by which the linear size L 
over which the periodic boundary conditions are imposed 
is taken to infinity or, equivalently, the lattice spacing is 
taken to zero, each first Chern number is quantized. 

The IQHE is an example in 2D for which the condition 
of constant Berry curvature F^^ ^ is met. In this context. 



the closed algebra (2.70a I was found by GMP (in fact, the 
algebra closes to all orders in q in this case) With the 
help of this algebra, GMP argue, within a single-mode 
approximation, that FQH states are incompressible. 

Recently, it was shown that lattice models with flat 
bands and nonzero Chern number also support incom- 
pressible FQH ground states,'^^^^ even though their 
Berry curvature is not constant over the BZ. This re- 
sult suggests to approximate the algebra (2.69a) by the 



algebra (2.70a), that is, to replace F^^, ^ with its average 
value over the BZ.SaHSI 



We will now consider a fully antisymmetric product of three projected density operators, expand it to third order 
in the momenta, and apply the same approximation by replacing each term in the expansion with its average over 
the BZ. The expansion reads 



^ f^Qi f^Qj f^Qk 



fceBZ I 



N 



\a=l 



I \^k+q,+q.,+q^ 



r 



(2.71) 



where the summation convention over the repeated in- 
dices i,j,k = 1,2,3 and J^, A = 1,2,3 is implied. The 
term of second order in q comes again multiplied by 
the Berry curvature, i.e., the density associated with 
the topological invariants Ch^ for A = 1,2,3 defined in 
Eq. (2.70b). As for the second term on the right-hand 



side, we recognize the integrand of the Abelian Chern- 
Simons form. We thus define 



^ E ^'^'"^P-.'^^A,fc- (2.72) 
fceBZ 



The sums 



N 



d^d^A^^k (2.73) 



fceBZ 



fceBZ 



vanish in the thermodynamic limit, since their integrands 
are derivatives of functions that are sufficiently well be- 
haved to allow for the application of the fundamental 
theorem of calculus. The term 



/i,fc 



fceBZ 



(2.74) 



for any triplet /i, J^, A — 1, 2, 3 is nonuniversal and in gen- 
eral nonvanishing. It vanishes if the system satisfies the 
isotropy condition 



(2.75) 



fceBZ 

for any triplet ijl,v,\ — 1, 2, 3. 



Finally, when replacing each factor in Eq. (2.71) by its average over the BZ, one obtains the approximation 

1 r 



Pq^ Pq, 



2ti 



i (g,: A g,) • Ch - (Qi A 92) • 93 - + ^"'q^c^q^ I^,, E 



"fc / \"fc+qi +92+^3 



fceBZ 



(2.76) 



where the summation convention over the repeated in- dices fc = 1,2,3 and fi,iy,X — 1,2,3 is implied. 
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This is the main result of Sec. Ill Cl Insulators for which 
the invariants Ch^ with A = 1,2,3 are nonvanishing can 
be viewed as a 3D extension of an IQHE or a layered 
system of 2D Chern insulators. In this case, Ch^ with 
A = 1,2,3 parametrizes the quantized off-diagonal part of 
the conductivity tensor.'^ The physics of such insulators 
is not intrinsically 3D and they are thus not our primary 
interest here. 

Even if the Berry curvature vanishes on average in the 
BZ so that Ch = 0, 6* can be nonzero and may take any 
real value in general. The value of 6 has measurable 
consequences as it contributes to the magneto-electric 
coupling in a 3D band insulator.'^ For 3D band insula- 
tors with either spin-orbit coupling that are time-reversal 
symmetric (symmetry class All) or with chiral symme- 
try (symmetry class AIII), 6 is restricted to integer mul- 
tiples of TT and represents a topological invariant.!^ We 
thus conclude from Eq. (2.761 that for 3D tight-binding 
Hamiltonians within the symmetry classes All or AIII, 
the fully antisymmetric product of three projected den- 
sity operators is dominated by the value of their topolog- 
ical invariant 9, just as the commutator of two projected 
density operators is dominated by the value of the Chern 
number in 2D tight-binding models within the symmetry 
class A. We will illustrate this statement with the help of 
a microscopic lattice model belonging to the symmetry 
class AIII in the following section. 



III. NONINTERACTING THREE-BAND 
TIGHT-BINDING MODEL 

The goal of this section is to define a "simple" single- 
particle Bloch Hamiltonian that supports a dispersion- 
less isolated band with nontrivial topological character, 
such that the electronic density operators projected to 
this ba nd ob ey the noncommutative algebra described 
by Eq. ( [2J6| with Ch^ = for A = 1, 2, 3 and 9 = tt. 
Our model belongs to symmetry class AIII and has three 
bands, which is the minimum number required to realize 
the desired 0-term.l^ One of the three bands is necessar- 
ily dispersionless as a consequence of chiral symmetry. 
Therefore, it can be taken as the basis for the construc- 
tion of fractional topological states in 3D. 



A. Definition 

We consider spinless electrons hopping between the 
sites = (r]^,r2,r3) of a 3D cubic lattice A and on- 
site orbitals, whereby each site r can accommodate three 
orbital degrees of freedom that we label with the Greek 
index a = 1, 2, 3. To accommodate the hybridization be- 
tween any of the three orbitals, we need to choose a basis 
for all 3 X 3 Hermitean matrices. We denote the unit 
3x3 matrix by Aq which, together with the eight trace- 
less Gell-Mann Hermitean matrices A„ with n = 1, • • • ,8, 
form the desired basis of all 3 x 3 Hermitean matrices. 



The second quantized tight-binding Hamiltonian is then 
defined by 



3+i 



-At) a 



M J2 Cl 



(3.1) 



r6A 



where we have introduced the 3-component operator — 
^c^.]^, cj..2, with 'cl a creating a spinless fermion at 
site r in the orbital a — 1,2,3 and obeying periodic 
boundary conditions under the translation r — >■ r -t- Lbj 
for any of the three orthonormal unit vectors e^, 63, and 
63 that span the cubic lattice A. This single-particle 
Hamiltonian depends on the real-valued parameter M . 

Translation invariance allows to diagonalize Hamilto- 
nian (3.1 1 upon performing a Fourier transformation on 



the creation and annihilation fermionic operators. If we 
denote with BZ the Brillouin zone of the 3D cubic lattice 
and with k any Bloch wave number from the BZ that is 
compatible with the periodic boundary conditions, then 



^=E^I^fe4 (3.2a) 

keBZ 

with the momentum-resolved single-particle 3x3 matrix 

(3.2b) 



4 



.7 = 1 



that depends on the 4-component real-valued row vector 

dk ^ {di^ l,d^ 2i'^k,3T'^kA) 

(3.2c) 



sin ki , sin k2 , sin fcg , M 



3 

E 



COS fc. 



With the help of the explicit representation of the eight 
Gell-Mann matrices from Appendix [B| one verifies that 

CHkC-^ ^-Hk, VfeeBZ, (3.3a) 

if and only if the 3x3 matrix C is given by 

C :=diag(l,l,-l). (3.3b) 

The fact that anticommutes with C implies that any 



pair of eigenstate u^^"* and u^. ' of Hj^ with nonvanish- 
ing eigenvalues are associated with the opposite single- 
particle eigenenergies e^^'' = — respectively. Since 
is a 3 X 3 Hermitean matrix for any wave number k 
from the BZ, it then follows that at least one eigenstate 



(-) 



must have the vanishing eigenvalue 



e^"' - 



(3.4a) 



irrespective of the Bloch wave number k in the BZ. For 
any Bloch wave number k in the BZ, the values taken by 
the nonvanishing eigenvalues 



-(+) 



-(-) 



(3.4b) 
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FIG. 2: Energy spectrum of the lattice model defined in 
Eq. (3.21 for different values of the parameter M. Panels b) 



and d) show the gap-closing topological transitions. Note the 
dispersionless band at zero energy in each spectrum. The 
spectrum is plotted along the straight path connecting the 
following points in the BZ: V = (0,0,0), X = (0,7r,0), M = 
(tt, tt, 0), and R — (tt, tt, tt). 



of M, i.e., whenever \dfj\ is nonvanishing over the entire 
BZ, there are two dispersive bands whose Bloch states 
are related by the chiral transformation 



,(+) - 



C 



,(-) 



and one dispersionless band u 



(0) _ 



C u"^^ of zero modes. 
Hamiltonian (3.2) breaks time-reversal symmetry, for 



the first three components of are odd wfiile the fourth 
component is even under k — )• —k for any value of M . 
This leaves no room for a particle-hole symmetry by 
which Hamiltonian (3.2) would anticommute with an an- 



tiunitary operator. Adding to Hamiltonian (3.2) any lin 



ear combination of the remaining Gell-Mann matrices X^, 
A2, A3, Ag, and the unit 3x3 matrix Ag breaks the chi- 
ral symmetry. Such perturbations change the symmetry 
class of Hamiltonian (3.2) from AIII to A. Although a 



chemical potential (a nonvanishing constant term pro- 
portional to the unit matrix Aq) does break the chiral 
symmetry, it does so by moving rigidly the entire en- 
ergy spectrum up or down in energy while leaving the 
Bloch states unchanged. The topological attributes of 
the three Bloch bands are thus untouched by the addi- 
tion of a chemical potential. 



follow immediately from the fact that any pair of the 
four Gell-Mann matrices A4, A5, Ag, and A7, anticom- 
mute pairwise while any one of these 4 Gell-Mann ma- 
trices square to either diag(l,0, 1) or diag(0, 1, 1). The 
minimum value reached by the magnitude |c?^| over the 
BZ thus determines the energy gap between the disper- 
sioness band of zero modes and the pair of bands related 
by the chiral transformation C. This energy gap depends 
parametrically on M and is nonvanishing if and only if 
\M\ 7^ 1,3. In turn, the corresponding Bloch states are 
derived as follows. One observes that the 2-component 
complex- valued row vector 



.1 



' (4a 



id 



k-A 



(3.5a) 



of unit length (q^ (7^ ~ 1) enters according to 




(3.5b) 



One then verifies that 



,(±) - 




,(0) _ 



Mi. = - 




(3.5c) 



are orthonormal Bloch states of 'H^ for any Bloch wave 
number k from the BZ. For any value of the parame ter M 
enter ing the single-particle Hamiltonian H, Eqs. (3.5c) 
and (3.4) define globally over the entire BZ the desired 
Bloch states with their dispersions. For generic values 



B. Topological invariants 

We shall take the thermodynamic limit L — >■ 00 with 
L the linear extend over which periodic boundary condi- 
tions are imposed. In this limit sums over wave numbers 
in the BZ are replaced by integrals over the BZ while the 
index k becomes the argument of functions. From now 
on, we shall identify the BZ with . We can then distin- 
guish two related topological invariants associated to the 
family of single particle 3x3 matrices T-L{k) labeled by 
the wave number k from a BZ with the topology of the 
3-torus owing to the periodic boundary conditions. 

The first topological attribute characterizes the bun- 
dle of Hamiltonians 'H(fc) over the BZ . For any wave 
number k ^ T^, there is a one-to-one correspondence be- 
tween the 3x3 Hermitean matrices H{k) and the vector 



d{k) e 



For any wave number k G T^, the mag- 



nitude |d(fc)| measures the momentum-resolved energy 
separation between the zero mode u^^^ (fc) and the lower 
and upper modes u^~^(fc) and u(^-'(fc), respectively. The 
eigenstates M'°^(fc), u^~\k), and u*^+-'(fc) are independent 
of the magnitude of |d(fc)|, i.e., they only depend on the 
coordinate defined by the unit 3- vector d{k)/\d{k) \ on the 
3-sphere . It follows that the topological attributes of 
the three Bloch bands of Hamiltonian (3.2) are deter- 



mined by the homotopy group Z of the map defined by 



keT-* 



dik)/\dik)\ e 



(3.6) 



between the BZ and the 3-sphere 5*^. For each pa- 



rameter M 7^ ±1,±3 entering in Hamiltonian (3.2 1, the 
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integer value taken by the topological invariant 



127r2 



T3 



(3.7a) 



determines which homotopy class the map (3.6) belongs 
to. Here, we are using the short-hand notation d^dj = 
ddj/dk^, with fi, X labeling the three coordinates of the 
momentum k and k, I labeling the four components of 
the vector field d, and the convention for summation over 
repeated indices. Explict computation of as a function 
of M delivers 



-2, |Af|<l, 
iy{M) = <j -1, 1< \M\ < 3, 
0, 3<|M|. 



(3.7b) 



Whenever = 1.3, the gap over the BZ closes at 
the discrete points (the lattice spacing is unity) 



TT (1,171,11), l,m,n^ 0,1. 



(3.8) 



These eight wave numbers change by a reciprocal wave 
vector under the operation of time reversal, under which 
k — > —k. In this sense, they are time-reversal invariant. 
The touching of the upper and lower dispersions at the 
wave numbers (3.8) occurs at zero energy and delivers a 
Dirac dispersion in their close vicinity when \M\ — 1,3. 



Hence, we call the wave numbers (3.8) Dirac points when 
\M\ = 1,3. For small deviations away from \M\ = 1,3, 
a spectral gap opens up at the wave numbers (3.8) that 
can be associated with a Dirac mass. Remarkably, the 
number of Dirac points that change the sign of their mass 
across a transition tuned by changing M through any 
one of the values |Af | = 1, 3 is equal to the change in the 
topological invariants (3.7). To see this, observe that the 



momentum resolved Dirac masses are given by 



du .A 
"'000'^ 


= Af - 3, 








.A 

"-001 


— du .A — 


di, .A 


= A/- 


1, 


dji, .A 




dv. .A 


= A/ + 


1, 


dr, .A 


= Af + 3. 









(3.9) 



With the help of these 8 integers, we define the integers 
'^b{M):=1 (-ir+"+'sign4_^^4. (3.10) 

m,n,/— 0,1 

The factor (_i)™+"+' assures that the mass sign is taken 
relative to the chirality of the kinetic piece of the Dirac 
operator. One verifies that (see also Appendix [d| 



(3.11) 



for any \M\ ^ 1,3. 

The second topological attribute characterizes the bun- 
dle of Bloch states u('^)(fc) over the BZ T'^ for any of the 
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FIG. 3: (Color online) Numerical evaluation of the topolog- 
ical invariant 6''-'^' = ■kv{M) (solid line) for a) the model (3.2) 
and b) the more isotropic model (3.21 1 with — 2. The 
anisotropy parameter I{M) defined in Eq. ( |3.18[ ) is equal to 
the difi^erence between the dashed and dotted line and van- 
ishes for specific values of M in panel b). The parameters 0' 
and 6" that sum up to the topological invariant 6 are defined 
in Eqs. (Glial and (Gild I, respectively. 



three bands d = — ,0, Whenever lAfj 1,3, it is 
nothing but the triplet of Berry phases^i 



e^^\M) ^ / d^ke^""^ Af') d^Af\ (3.12a) 



T3 



where we have introduced the Abelian Berry connection 

4'^)(fc):=(^u('^)tA^(")^(fe) (3.12b) 

for any of the three bands a = — , 0, +. With the help of 
Eq. (3.5c), one deduces that 



0(o)(Af) d^fce^"^ {q^ d^qd.q^ d^q) (fc) (3.13a) 
and 

0(-)(A//) = 6l(+)(Af) = ^e'-°\M) (3.13b) 

when \M\ ^ 1,3. Explicit evaluations of the Berry phase 
of any of the three bands then yields 

e'^-\M) = e^+\M) = ^e^°\M) = ^i^{m) (3.14) 

when \M\ ^ 1,3. 



C. Density algebra 

We will now argue that the algebra obeyed by the 
momentum-resolved density operators 



fceBZ a=l 



(3.15) 
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after projection to the isolated flat band of Hamilto- 



nian (3.2 1 is dominated by the integrand of the topo- 
invariant = tt in the long wave-length and 
For that, we use Eq. (2.761 and note 



logica 
low energy limit. 



2.70b) are 



that the Chern numbers Ch defined in Eq 
vanishing for all bands of the model defined by Eq. (|3.l|) 
(see Appendix [Cj 



The values taken by '^^^x defined in 
For the g auge 
we 



Eq. (2.74) are in general gauge dependent. 



choice of the Bloch eigenstates used in Eq 
show in Appendix [C] that 



(3.5c 



case of a partially filled Landau level. Henceforth, one 
may expect interesting many-body ground states to ap- 
pear once electron-electron interactions are added to the 
model. In that regard, we observe that any many-body 
Hamiltonian that includes an interaction build out of the 



projected density operators (3.15) is invariant under the 



chiral transformation (3.3a), since the projected density 



operators (3.15) themselves are invariant under the chiral 



transformation (3.3a 



0, 



/^,i/ = l,2,3, (3.16) 



D. 



Surface states 



and 



while 



^312 



^321 



^132 ^ ^12 



^213 



-I{M)/2, 



(3.17) 



(3.18) 



Here, I{M) is some nonuniversal real number, that is 
nonvanishing generically and measures by how much 



(3.19) 



fceBZ 



fails to be proportional to e^^.^- Equation (2.761 thus 
reads for the fiat band of the 3-orbital model 



JJk-^o) -^0) -^0) 
f^Qi f^Qj ^Qk 



^1 
feeBZ 



,(0) 



fe+qi+q2+'?3 



(3.20) 



It turns out that I{M) is nonvanishing for any nonvanish- 
ing value of M in the model defined by Hamiltonian ( 3.2 1 



This can, however, be remedied by replacing the vector 
in Eq. (|3.2c|) with 



: (sin ki , sin , t] sin , 
M — cos ki — cos k2 — f] cos fcg). 



(3.21) 



where the real parameter 77 was introduced. While rj = 1 
corresponds to the original model, for large enough 77, 
one can find an AI = M{r}), such that 



/(77,M(r7)) =0. 



(3.22) 



In Fig. |3j the case 77 = 2 is exemplified. 

In summary, we have found a 3D lattice model with 
a topologically nontrivial fiat band. We have shown 
that the electron density operators, when projected to 
this flat band, obey a noncommutative algebra that is 
dominated by the value of the topological index of the 
model. Upon partial filling, this flat band thus provides 
a manifold of many-body noninteracting ground states 
with macroscopic ground state degeneracy, similar to the 



We shall here provide an interpretation of the topo- 
logical invariant (3.10) as a manifestation of the surface 



states associated with a spatial dependent mass parame- 



ter M in the Hamiltonian (3.2). This observation applies 



when considering the surface states that connect bands 
separated by a bulk gap. Such surface states, connecting 
the upper and lower band, appear only when the periodic 
boundary conditions are replaced by open boundary con- 
ditions that implement a slab geometry with the surface 
normal parallel to the direction. 

In order to study the surface modes, we consider the 



low energy description of the Hamiltonian (3.2) by lin- 
earizing it around each of the 8 nodal points in the Bril- 
louin zone fc^ ^,„ — n{l,Tn,n), with l,m,n = 0,1. The 



Hamiltonian (3.2 1 in the linearized approximation fac- 



torizes according to 



n= (g) % 



(3.23) 



l.m,n— 0.1 

For example, the expansion around fcooo produces 



K 



000 



k_ 

fcg-iA/oool, (3-24) 



i fc I fco + i AI, 



000 







where kj. = ki± ifcj, fc - 



-ia,, ,, for j 



1, 



, 3 and 



000 = -^'^ ^ 3. For a uniform mass Mqqq, the spectrum 



breaks into three low energy bands with eigenvalues 
and ±Vfc2 + |Mooo|2. 

We now regard Mqqq as a domain wall conflguration 
along the rg-direction, which we choose to parametrize 
as 



(3.25) 



where 8 is the Heavi side function. The choice of a sharp 
domain wall in (3.25) facilitates the analytic treatment of 



the eigenmode equations and does not affect the general- 
ity of the following discussion. Due to the translational 
invariance on the ej^-ej plane, we seek solutions of 



"^000 ''Pi 



000 V^OOO.K 



000, K I 



(3.26) 



with %oo J p^r^) 



"^000 



^(rg), whereby p 



{ri,r2) and k = (ki,k2) are, respectively, the coordi- 
nates and momenta projected on the plane. The 
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components of the spinor wavefuction 

</'000,k(''3) (/«(^3)>5«(^3):^«('"3))^ (3-27) 

satisfy 



(3.28a) 



iMooo(^3) h^ir^), (3-28b) 



and 



At ^ 0, the solution of Eq. (3.28c) yields 



(3.28c) 



(3.29) 



where is a normalization constant and A ^ := 
Mqqq + /«2 — £^ > 0, while the delta function discon- 
tinuity at r3 = imposes the condition = Mqqq. 



Therefore, the domain wall configuration (3.251 bounds 
surface states with dispersion 



(3.30) 



provided Mqoo > 0. Evaluating the solution (3.29) 



in (3.28a) and (3.28b) yields the spinor wavefunction, 



which, up to a normalization constant J\f, reads 



(3.31a) 



'i^OOO,±,K; 



,0,1) 



K 



(3.31b) 

The discussion of the boundary states of the low energy 
Hamiltonians with n — 0, Hi^q, is very similar to that of 
7^000- ^^^^ case, the exis tence of gapless surface states 
with dispersion as in Eq. (3.30) for sharp domain wall 
configurations 



(3.32) 



requires M > 0. The explicit form of the eigenspinors 
(omitting the dependent part) is 



^ioo,±,. = 2-^/^ (Te+'"^0,l)^ 
^m±.. = 2-^/' (Te-'"^0,l)^ 
Voio,±,. = 2-'/' (±e+'"^0,l)^ 



(3.33) 



For the boundary states of the low energy Hamiltoni- 
ans with n — 1, T-Liyjii, the extra minus sign coming from 



the Taylor expansion around ^3 = tt implies that the gap- 
less surface states exist for domain wall configurations 



(3.34) 



provided Af;„ii < 0. The eigenspinors in this case read 

¥'ooi,±,. = 2-1/2 (±e--»,0,l)\ 
^ioi,±.. = 2-1/2 (^,+io„^o,l)\ 
'/'oii.i,. = 2-1/2 (±e+-»,0,l)\ 
'/'iii.±.. = 2-1/2 (Te--»,0,1)\ (3 35) 

In order to account for all the possible surface modes 
in a finite size configuration, we now take, for the sake 
of concreteness, our system to be a slab, infinite in the 
e]^-e2 plane and confined in the rg-direction by 



< 



< r^°"°™, with r3°P - r5°"°"^ assumed to be much 



„bottom 
'3 



larger than any other length scale so as to regard the 
two surfaces as completely decoupled from each other. 
Moreover, let us adopt the convention that the vacuum 
is characterized by a positive value of the gap parameter 
(Afyjjj. > 0), which then changes to negative values for 
r^^^ < < T-bottom -poY this particular configuration, 
the discussion above implies the presence of gapless sur- 
face states associated with "H^^q {'M-imi) at the surface 

^3 = rf^ {r3 = r^°''°n for M^^o < (M,„,i < 0). 

In order to make a connection with the topological 
invariant (3.10) we now compute the winding number of 



the eigenspinors as 



1^1, 



ImO 



d/t • ((y3|,„o_^_^V^ Vi,no,±,n) : (3-36a) 

) , (3.36b) 



where the explicit overall sign difference between (3.36a 



and (3.36b) reflects the opposite orientation of the out- 
ward normal vectors -\-e.^ and 



and — r- 



bottom 
3 



-63 at the surfaces 



^, respectively. Direct computation 
using Eqs. (|3.33[) and (3.35) gives 



'100 



-1, 
-1. 



(3.37a) 
(3.37b) 



The total winding number of the surface states is en- 
coded in the quantity 



E 



(3.38) 



M,, 



, <0 



which acquires the following values: 



V = < 



+2, \M\ < 1, 
-1, 1 < \M\ < 3, 
0, 3 < \M\. 



(3.39) 
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Comparison between Eq. (3.39) and Eq. ( 3.7b[ ) thus es- 
tablishes a direct relationship between the topological in- 
dex (3.101 and the total winding number of the surface 
states i>. Similar analysis of the finite size system spec- 
trum for domain wall configurations of the gap parameter 
along either the x or the y directions reveals the inexis- 
tence of surface states. 



IV. INTERACTIONS WITHIN THE 
SINGLE-MODE APPROXIMATION 

We begin by reviewing the single-mode approximation 
(SMA) to the FQHE from Ref. EH 

In the IQHE, the external magnetic field organizes 
the single-particle spectrum into degenerate Landau lev- 
els, whereby two consecutive Landau levels are sepa- 
rated by the energy gap Hto^. The cyclotron frequency 
oj^ = h/lm^l'^g) is proportional to the magnitude B of 
the uniform magnetic field. 

We consider the limit of very strong magnetic fields rel- 
ative to the characteristic energy scale V of the electron- 
electron interactions, i.e., hiu^ 3> V. Moreover, we con- 
sider a filling fraction = ^/^q < 1 (<i> the magnetic 
flux and $q the flux quantum) such that the exact many- 
body ground state j^'g) does not break spontaneously any 
symmetry. The translation invariant interacting Hamil- 
tonian H describing a nonvanishing density of spinless 
fermions moving in a plane perpendicular to an external 
magnetic field of uniform magnitude B and interacting 
pairwise with a (screened) Coulomb interaction is then 
well approximated, as far as low energy properties go, by 
its projection H^^^ onto the vector space spanned by the 
lowest Landau single-particle levels. 

Upon imposing periodic boundary conditions in an 
area of linear size L, ^lll given by 



where 



(4.1a) 



(4.1b) 



is the Fourier transform of the screened Coulomb inter- 
action, while 



'5pq:=Pq-(«'o|pJ*o) 



(4.1c) 



is the Fourier component of the fermion density operator 
after projection into the LLL measured relative to its 
expectation value in the exact many-body ground state 

l*o)- , 

Inspired by the early work of Feynman and Bijl in their 
study of excitations in ''He,'^^ CMP in Ref. [21] consider 
the variational state 



(4.2) 



whose energy expectation value A^, measured relative 
to the exact ground state energy Eq, sets a variational 
upper bound on the low exc itation spectrum of the LLL- 
projected Hamiltonian (4.1). 



Assuming the inversion symmetry 



(4.3a) 



a direct calculation using Eqs. (4.1) and (4.2) leads to 
Afc = ^, (4.3b) 

where 



/fe 



5p_ 



LLL J 



Sp. 



+k 



and 



(*o|'5p-fc^P+fc|4'o) 



(4.3c) 



(4.3d) 



One recognizes on the right-hand side of Eq. (4.3d) the 
static structure factor. The insight of GMP in Ref. i3T] 
was to realize that the density operators projected onto 
the lowest Landau level close the exact algebra 

[Pq, Pk] = 2i sin ( (q- X fe) • 63 4 ) pg+fc (4.4) 



{£g is the magnetic length). In turn, the algebra (4.4) 
implies that 



fk - \k\ 



(4.5) 



in the small |fc| limit. Hence, in the FQHE, a necessary 
(but not sufficient) condition for the existence of a finite 
gap in the thermodynamic limit is to have 



Ifcl' 



(4.6) 



also hold in the small |fc| limit. In fact, Eq. (4.6) was 



shown in Ref. [311 to be satisfied when \'^q) is chosen to be 
any Laughlin state with filling fraction h' — 1/m, where 
m is an odd integer. 

In the spirit of GMP, our starting point is a single- 
particle Hamiltonian defined on a c?-dimensional Bravais 
lattice and sharing its point group symmetry. We also 
assume that there exists at least one band that is inde- 
pendent of the lattice momentum, i.e., a flat band, and, 
furthermore, that is separated from the other bands by 
a single-particle gap A. We constructed a 3D example 
thereof in Sec. [III| We then imagine switching on adia- 
batically a pairwise interaction that preserves the Bravais 
lattice point-group symmetry, say a (screened) Coulomb 
interaction. We shall denote with V the corresponding 
characteristic interaction energy scale. In the regime for 
which V, Hamiltonian (4.1) can be reinterpreted as 



the interacting Hamiltonian projected onto this flat band, 
provided we identify as the Fourier transform with the 
lattice momentum q of the pairwise fermion interaction. 
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as the Fourier transform with lattice momentum q 
of the projected operator describing density fluctuation 
measured relative to the fermion density with lattice mo- 
mentum q of the exact many-body ground state |5'q), 
whereby we assume that j^'p) does not break sponta- 
neously any point-group symmetry of the lattice. 
The projected density operator on a flat band reads 

^ _ \ " t - 

Pk — 2^ ■ "p+fc Xp Xp+k 



^'^p,k Xp Xp+fc) 



(4.7) 



where Uj^. g C is vector-valued (its components range 
over the number N of orbitals per site of the Bravais lat- 
tice), while Xk '"^^^ Xk the annihilation and creation 
operators, respectively, of single-particle fermionic eigen- 
states on the isolated flat band with lattice momentum 
k. Hence, they satisfy the canonical fermionic anticom- 
mutation relations 

{xfc^Xfc'} = {xLxL} = 0, [xk,xi'} ^ S^j,, (4.8) 

for any pair k and fe' of lattice momentum. In carrying 
out the program laid out in Eq. (4.3) for a general lattice 
Hamiltonian with a flat band, one notices two immediate 
obstacles. 

The flrst one arises from the fact that the commutator 
of two ( proje cted) density operators does not satisfy the 
algebra (|4.4|) found by GMP for the FQHE in a uniform 



magnetic field. However, it was noticed in Ref. ,49. that. in 
the limit of small lattice momenta k and k' , the commu- 
tation relation between two projected density operators 
reads 



[p(fc),p(fc')] = / i {kAk')-B{p 



(4.9a) 



X XHP) x{p + k + k') 



in the thermodynamic limit i — > oo, whereby the short- 
hand notation 



{2-k/LY 



(4.9b) 



is the (real-valued) Berry field strength of the fiat band, 
and 



A{p) 



V u) (p) 



(4.9d) 



is the (imaginary-valued) Berry connection of the flat 
band, while • • • in Eq. (4.9a I accounts for higher order 
terms in powers of k and k' . Consequently, it was pro- 
posed in Ref. 132! that the numerical observation of the 
FQHE without an external magnetic field in 2D Chern 
insulators in Refs. I45H481 can be understood on the ac- 
count that, because in a 2D Chern band insulator the 
integral of the Berry curvature on the Brillouin zone 
equals the (nonzero) Chern number, replacing t3{p) in 



Eq. (4.9 1 by its average, implies the GMP algebra (4.4) 



in the long wavelength limit. However, we would like to 
stress that, contrary to the 2D Chern band insulators for 
which one can associate the notion of an average Berry 
curvature due to the nonzero Chern number, for the 3D 
lattice models studied in Sees. |lT] and the integral of 
the Berry curvature vanishes so that replacing 3{p) by 
its average is meaningless. Even for 2D Chern band in- 
sulators, the Berry curvature is generically nonuniform; 
a fact that should be reflected in the exact many-body 
wavcfunction. 

The second obstacle to applying the SMA to an inter- 
acting lattice model is the fact that no good candidate 
wavefunction is presently known with which one can com- 
pute the static structure factor Sj^ and compare its small 
k dependence with that of as was done by GMP in 
Ref. |3T1 Nevertheless, information about the behavior 
of for small k and the requirement of a finite gap in 
the thermodynamic limit, i.e., Aq 7^ for fe — > 0, 

puts a constraint on the static structure factor for small 
fe and, correspondingly, on the correlations of the exact 
many-body wavefunction. 

In Appendix IE] we discuss in detail the evaluation of 



is used, 



the function f^. defined in Eq. (4.3c I to lowest order in 
fe. Our main result is that, due to the nonclosure of the 
density algebra for any d-dimensional lattice model, the 
leading contribution to /(fe) reads 



3{p) := -i (V A A) [p) 



(4.9c) 



/(fe) = j J J v{q)[{kAq)-5Bip)\[{kAq)-SB{p')\^h{p)nip'] 

q p p' 

v{q) ^ (fe A q) • (d^B) {p)k^^ [( 5p{-q) xHp) x{p + ?)) - ( xHp + Q) x{p) Sp{q) 



q p 



(4.10a) 
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where J = J d'^q/{2-K/LY and the summation conven- 
q 

tion is imphed over the repeated index /i = 1, • • • , rf. In 
Eq. (lilOal), 



while 



n{p):= x'{p)x{p) 



(4.10b) 



with chiral symmetry. We estabhshed a relation between 
the algebra of the projected position operators and the 
classical Nambu-bracket of volume-preserving diffeomor- 
phisms of 3D fluids, that might bridge the description of 
classical ideal fluids and that of topological incompress- 
ible states in 3D. 



SB{p) := 3{p) - 3 



(4.10c) 



denotes the deviations of the Berry curvature 3{p) away 
from the uniform background value 3. This uniform 
background value is defined in such a way that, when 
d = 3, 



Ch^ :=27r X 



d^p 



2 J (27r)3 
Z3 ^ 2 y (2^ 

T3 



B\p) 



27r 



1 



(4.11) 



B 



with A = 1,2, 3 is compatible with a generalization of 
the 2D Chern number to layered (quasi-2D) materials. 
The result (4.10a) should be contrasted with the calcu- 
lation in Ref. |31j for which the order term in /(fe) 
vanishes identically as a consequence of the algebra (4.4 1. 



The formula (4. 10a I thus establishes a direct relationship. 



within the SMA, between the deviations of the Berry field 
strength away from a uniform configuration and the order 
fe^ contribution to f{k) 



f{k)^\k\ 



(4.12) 



Such a relation is relevant either for 2D fractional Chern 
band insulators for which, despite a nonzero Chern num- 
ber, 3{p) can be nonuniform throughout the Brillouin 
zone or for the general classes of 3D lattice models stud- 
ied in Sees, [n] and III for which the integral of 3{p) van- 
ishes. The result (4. 10a I also indicates that a prerequisite 



for the existence of a nonvanishing but finite many-body 
gap to excitations above the many-body ground state 
is that the static structure factor s{k) has also to van- 
ish as fc^ to allow for the possibility of a nonzero ratio 
A(fc) = f{k)/s{k) and therefore a nonvanishing SMA 
gap in Eq. (|4.3b). 



V. SUMMARY 

The noncommutativity of coordinates and density op- 
erators in a featureless liquid-like electronic state can be 
a local probe of its topological character. In this paper, 
we have studied how this fact, which is well-established 
for quantum Hall fluids in 2D, carries over to 3D topo- 
logical states of itinerant electrons. In the limit of long 
wavelength, we found that both the algebras obeyed by 
projected position and density operators are character- 
ized by the topological invariant of a 3D band structure 
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Appendix A: Chern-Simons invariant as product of 
covariant derivatives 

In this Appendix, we will review a derivation sketched 
in Ref. 59, that relates the trace over the 3- bracket of 
covariant derivatives to the Chern-Simons invariant in 
3D. We shall use a representation in terms of differential 
forms, which is equivalent to the notation using the fully 
antisymmetric 3- bracket used in the main text. Define 
the covariant derivative to be 



D^:=a, + A„ M=l,2,3, 



(Al) 



where A is some (anti-Hermitean) non-Abelian gauge 



field. Define the one form 



D := D„ dx^ 



(A2) 



and the two form 

F := (d^ + A^) [d^ + A^) dx'^ A dx" = DD. (A3) 
We begin with the manipulation 
F =DD 

= (a^ + Aj {d, + A,)dx^^ ^dx^ 
= (d^d, + d^A, + A^d, + A^A^i dx^ A dx'' 



= + {d^A^) + A^d^ + A^d, + A A dx^ A dx"" 



( 



{d^A,) -A^d, + A^d, + A^A,) dx^ A dx^ 
[{d^A,) + A^A,] dx^^^dx'' 



(A4) 



where it was assumed that the derivatives are commuta- 
tive 

{d^d, - d,d^) ■ = (A5) 
when acting on everything to the right. 
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The 2-fomi F has no derivatives acting to the right 
anymore, i.e., it is not a differential operator when acting 
on the basis that diagonalizes the position operator. We 
thus know how to compute the trace TrF of F. 

The situation is different for the 3-from arising from 
the third power of the covariant derivative 



= DF. 



(A6) 



In general, it has derivatives acting on the right. In order 
to take the trace of this quantity, we are going to rewrite 
these derivatives as commutators of a 2-form with a 1- 
form. Indeed, upon taking the matrix trace of the com- 



mutator of a 2-form B and a 1-form C, 

Tr[B, C] - Tr {B^^C^ - C^B^^) dx^ A Ax^ A dx^ 

= Tr {B^,C^ - B,^C^) dx^ A dx" A dx^ ^^^^ 
= Tr {B^,,C^ - B^^C^) dx^ A dx" A dx^ 
= 0, 



where we have used the cyclicity of the trace to reach the 
second equality and made a cyclic permutation of the 
indices on the second term to reach the third equality. 



We now do the expansion 

DF = {d^ + A^) [{d^A^) + A^A^] dx^' A dx" A dx^ 

= [MdM + d^A^A^ + A^idM + A^A^A^] dx^ A dx'' A dx' 

Observe that the first term on the right-hand side of Eq. (|A8|) can be rewritten as 



(A8) 



d^{d^A^)dx^' A dx" A dx^ = {{d^d^As^) + a^^a^} dx" A dx" A dx^ 
u„„ = {{d^d^A^) + d^A.d^} dx^ A dx- A dx^ 

= { (d^d^A^) + A.d^d^ + [d^, A,d^] } d.;^ A dx'' A dx^ 
ti„„ = { (d^d^A^) + A^d^d, + [d^, A,d^] } dxf^ A dx" A dx^ 

= {d^AAx + [9^, Ad^] } d.T^ A dx-^ A dx^ 
^[d^,A^d^]dxi' Adx" Adx^, 

where we again used that d^d^^ — d^d^ = to set the underlined term to zero, while the second term on the right-hand 
side of Eq. ( A8 ) can be rewritten as 



cyclic pcrmutatii 



cyclic pcrmutatii 



(A9) 



1 



1 



1 



d.A^A^dx'' A dx'' A dx^ = - [a, , A^Ay] + - [d^ A^, Ay] + -A^A^d^ + -A^^d.A^ \ dx^' A dx"^ A dx 



= \ \[d^..AM + ^[9^^., A] + l^^dM \ da:^ A dx^ A dx^ 



(AlO) 



Thus, we have rewritten Eq. (AS) with all derivatives acting on the trace now represented by commutators between 
one and two forms, 



DF 



[9^, Kdx] + \[d^. A,A^] + \[d^A,,A^] + \A^,{dM + + A^A.A^^ dx'^ A dx'' A dx^ 



(All) 



According to Eq. ( A7 ) , we can discard all commutators when taking the trace of D'' 



Tr DF = Tr D'' 



- Tr ( -A^{dM + A^A.A^ ] dx'^ A dx'^ A dx^ 



(A12) 



Appendix B: Gell-Mann matrices 

The Gell-Mann matrices are 3x3 Hermitian matrices 
that are a representation of generators of SU(3). They 



are defined as 





(Bl) 
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and 




(B2) 



Appendix C: Density algebra for the three-orbital 
model 



In this Appendix, we evaluate the fuUy antisymmet- 
ric triple-product of projected density operators given by 
Eq. (2.761 for the projection on the dispersionless mid- 



dle band of the three-orbital model defined by Eq. (3.1) 
Specifically, we will compute the Chern numbers 



Ch-" 



(2^)2 



^y'd3feF^,(fc)eZ, (CI) 

BZ 



for A = 1, 2, 3, the Chern-Simons-invariant 



(C2) 



BZ 



and the term 



4(27r) 



^ J d^k [A,ik)d^ + A^{k)d,]A^{k) 



BZ 



(C3) 



for fijiy^X — 1,2,3 in the thermodynamic limit. (We 
have dropped the symbol refering to the projection for 
notational simplicity.) While both Ch and 6 are gauge- 
invariant, I^^A with iJ.,h',X — 1,2,3 depends in general 
on the gauge chosen for the Bloch basis of the single- 
particle Hilbert space. For the three-orbital model de- 



fined in Eq. (3.1 ), th e block off-diagonal projector q{k) 
defined in Eq. (3.5a) delivers a natural choice of gauge 



for the Berry connection of the flat band 

A{k) ^q\k)Vq{k). 
In this case, A can be decomposed as 



(C4a) 



A'{ki, ^2, fcg) -|- A"{ki, fca) 
Aik)= \ A'{k2,k„k,)-A"{k2,k„k,) I , (C4b) 

^3(^1' ^2; ^3) 



where 



A'{k) = -i 
A'\k) = +i 



sin fcj sin fcg 
G{k) ' 

cos ki sin ^2 
G(fc) ' 



(C4c) 
(C4d) 



_ . 1 + cos fcg (cos fci + cos fc2 - M) 
Asik) - -1 ,(L4e) 



G{k) = 3 + [ M - ^ cos fc^ j - cos^ k^ 



(C4f) 



It follows that 

A [ki, ^21 ^3) 

as well as 

A (ki, fc2, fcg) 



— A (— fci, k2, k^) 
+ A' (k-^, "^2' ^3) 

— A (fc]^, k2, ~k^), 



+ A ( — /C]^, fc2, fcg) 

— A [ki, —k2, fcg) 
+ A" [ki, fc2, — fcg). 



(C5a) 



(C5b) 



while A^{k) is an even function of ki, fcj, and k^. 

As a consequence, all terms appearing in Fi^{k) and 
-F23(fe) are an odd function of either k-^ or ^2. Thus, 



Ch^ = Ch^ = 0. 



(C6) 



Furthermore, 



Ch^ cx / d?kF, 
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BZ 



= j d'k{d^ [A'{k2,k„k,)^A"{k2,k„k,)] 
BZ 

- d2[A'{k„k2,k,) + A"{k„k2,k,)]} (C7) 

J d3fca2^"(fc2,fci,fc3) 
BZ 

^ ^2[A"i27T,k„k,)-A"{0,k„k,)] 
= 0, 

since 92^"(fc2, fc^, fcg) is a continuous function of ^2 with 
periodicity 27r. We conclude that the Chern numbers Ch 
defined in Eq. (CI) vanish identically. 



To calculate 6 andI^i,A, defined in Eqs. (C2 ) and ( C3 ), 
we consider integrals of the form 



J d'kA^d^A^, 



fi,iy,X^ I,- ■ ■ ,3. 



(C8) 



BZ 



By partial integration, one finds that these terms van- 
ish identically ii ^ — X. Using the symmetry relations 
Eq. ([C5| it follows that 



^M^3^3' ^3^3^, 



M= 1,2, 



(C9a) 



can be decomposed into terms that are either odd in k^ 
A^^d^A^, A^d^A^, M=l,2, (C9b) 



or odd in k^ while 
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can be decomposed into terms that are either odd in 
e'^'^^k^ or odd in and 

A.d^A^, A^d^A,, M = l,2, (C9c) 

are both odd in k^. Thus, the A; integrals over aU these 
quantities vanish. It then foUows that vanishes if at 
least any two of fj,, v and A are equal. 
It remains to consider the terms 

Jd'kA^d^A^, ti^v^K (CIO) 

BZ 

which are nonvanishing in general. On one hand, defining 



BZ 



StT j ^"(^1' ^21 ^3)^2^3(^1' ^2' ^3)' 

BZ 



(Clla) 



partial integration delivers 

-^' = hj ^'^^3^2^, (Cllb) 



BZ 



and using the identity A^{ki, fcj, k^) — A^{k2, k^, k^) one 
obtains 



1 

8^ 



d^kA^d^A^, 



BZ 



9' = ^! d'kA,d,A,. 



(Clfc) 



BZ 



On the other hand, defining 

+ 0" •= ^/ d'fe^^3^ 

BZ 

partial integration delivers 

-^" = ^/ d'kA,d,A,. 



BZ 



Finally, numerical evaluation of 

0^46' + 29" 



(Clld) 



(Clle) 



(C12) 



reveals that 9 is quantized in units of tt as announced, 
while 9' and 9" are not quantized and are not equal in 
general (see Fig. [s]). Furthermore, we find 
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-132 



2^231 — 21^213 — ^" 
2^312 — ^321 — 0, 

which after the identification 



/= -{9" ~9') 

TT 



(C13a) 

(Cf3b) 
(CI3c) 

(C13d) 



is the announced result. 



Appendix D: Equivalence of Chern-Simons and 
Dirac invariants 



The purpose of this Appendix is to prove that the 
Abelian Chern-Simons invariant, defined by 



1 

An 



d'ke^'^'A^d^A^, 



(DI) 



T3 



with the Abelian Berry connection A^(k) is equivalent 
to the Dirac invariant defined in Eq. (3.101 for the 
case of a Bloch Hamiltonian with chiral symmetry and 
three bands. The topological attributes of such a Hamil- 
tonian are characterized by its normalized off-diagonal 
part q{k) from Eq. (3.5a) in terms of which the Abelian 
Berry connection reads 



AJk)=q-\k)d^qik). 



(D2) 



Here, q{k) represents a map from T3 (the BZ) to and 
9/'K is the associated winding number. As a member of 
q{k) can be parametrized by three angular coordi- 
nates 



cos a e 



(. 

\ sm a e 



and the Berry connection reads accordingly 
A — i cos^ a d^ip -f i sin^ a d^-d, 



(D3) 



(D4) 



where we suppress the variable k for the moment. As we 
shall see, contributions to the winding number (Dl ) arise 
from vortex lines in <p(fe) and ??(&). Rewriting 



An 



d^k sin 2a (d^a) (9^i?) d^ip 



T3 



An 



d^k cos^ a {d^d) dyip 



d'k cos^ a [{d^d,^) d^^ + (d^d^p) a,^] 



T3 



(D5) 



the antisymmetric double derivatives in the last term 
contribute a delta-function for k on the vortex lines times 
the winding of the vortex. 



Let us now specialize on the model given by Eq. (3.1 ) 
in which case 



(fi ~ arg (sin k^ + i sin fcj) , 



arg 



sin k^ + i I M — cos k^ 



4=1 



(D6) 
(D7) 



and cos a = 1 in the vortex lines of ip, while cos a = 
in the vortex lines of Observe also that the first 
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term in Eq. (D5 1 vanishes, since the either of the partial 
derivatives d^'d and d^ip vanishes on each surface with the 
normal k^. The four vortex lines of f are parametrized 
by 



{nnT,mr,k^), m, n€E{0, 1} (D8) 



and their winding numbers are (—1)™+" . Eq. ( |D5[ ) then 
simplifies to 

1 ^ r 



1 ^ (-l)™+"+'sign4^„,^. 



m,n,l—0 



(D9) 



where we have written the number of phase windings of 
(p in the vortex line of t9 as 



dfcg 

2tt 



sigi^^fc„.„o;4-signd^,_^^^.4 



sign 4^ 



i=0 



(DIO) 



and is defined as in Eq. (3.8). In writin g Eq. (D9|, 
we have recovered the Dirac invariant (3.101. 



f2,k ■=l^^q( [^P-q^^P+k] [SP-k^Sp+q] )> (Elc) 



/3,fc ^ X! ( [^P-k^ [^P+q' ^P+k]] ), (Eld) 
q 

and 

Uk •= ^ H ( [^P-k^ [^P-q^^P+k]] ^P+q)- (Ele) 



The commutator of two projected density operators 

as 



can be expressed, with the aid of Eq. ( 4 



[Pq;Pfc] — -^p.q.fc Xp Xp+q+fc; (E2a) 



where 



^p.q^k ' ^^p.q^^p+q.k -^^p+k.q ^^p.k' 



(E2b) 



The nested commutators of three projected density op- 



erators can be expressed, with the aid of Eq. (4.8 1, as 



Pfc> Pk]\=^ ^p,qM Xp Xp+q, (E3a) 



Appendix E: SMA for a flat band 

We present some of the intermediate steps needed to 



derive Eq. (4.10a). (For ease of presentation, we use 
Latin instead of Greek indices for the momentum com- 
ponents in what follows. Summation convention over re- 
peated indices is also implied.) 



Our aim is to evaluate Eq. (4.3c) up to order q k . The 



commutator in Eq. (4.3c ) can be conveniently broken into 
four contributions, 

fk = fl,k + /2,fe + f^M + h,k^ (Ela) 
each of which read 

/l,fc ■=\^^q{ [^P-k^ ^P-q] [^P+q^ ^P+k] )> (Elb) 



where 

^p,qM ■= Rp-k,q,k ^p-k ~ ^p.q,k ^^p+q+k.-k' (E3b) 

Observe here that the identity 







t 




The 


Pfc> [Pq^Pk] 




P-k 



P-k, [P-q,P-k\ 



implies that 



A* = A 

p,q.k p+q. — q, — k' 



(E4) 



(E5) 



Needed is the expansion of ^ and up to 

order q^k^ . We start with 
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^'^p,q ^ ' ^p+q 



ip + q' a, Up+- q'q^ d, djU^ + 



1 + q'ul- d^Up + - q'q^ u], • d^d^Up H 

1 + q' A^^ p + ^ g'gJ' u], • 9, ^^Up + • • • 



(E6a) 



where we have introduced the (imaginary-valued) Berry connection 



^i,p — ^p ' ^p 



(E6b) 



and the summation convention over repeated indices i, j = 1, • • • , c? is imphed. The symbol with i = 1, • • • , d is to 
be regarded as a derivative with respect to the argument of the function on which it acts. Similarly, 



+ I {q'q'k"^ + k'k^q'^') ' d,d^d„,u + \ {k'k^q^^ + 2g"g'F) d„,u^ ■ d,d^d„,u + ^ q^q^k"^ d^d^u^ ■ d„ 
+ J q'q^k^k"' ■ a, djdi d„,U + \ k^yq^q"' 9; v) ■ O^d^u + ^ q'q^k^k"" 9,- d^v) ■ 5, d^d^u +■■■ 



(E7) 



where the summation convention over the repeated indices i,j,l,'m = l,--- ,d is implied. 

We multiply Eq. (E6al by Eq. (E7) and antisymmetrize with respect to the interchange of q and k. We obtain 



Rip, q, k) = (7;f ) (p) + (FFg™ - q^q 



.(3) 



(E8a) 



where the summation convention over the repeated indices i, j,l,m — 1, ■ ■ ■ ,d is implied and we have introduced the 
short-hand notation 



) ip) := {F,^) ip) ^ {d,A^ - d,A^) (p), 



and 



di d„^v) ■ d, u~d, dj • 9, a,„u + 2^, 9„ {u^ ■ d, u) - 2A^ dj (u^ • 9, d^u) 
+ 29; • di d^u- 2d, v) ■ 9; d^u (p) 
for i,j,l,m = l,---,d. We evaluate 

A(p, q, k) = R{p - fe, q, k) M (p, -fe) - i?(p, q, k)M{p + q + k, -k) 



(E8b) 



(E8c) 



(E8d) 



l-fc%(p) 



- Rip,q,k)-k-d^R{p,q,k) + - 

- i?(p, q, fc) - fc'^<9,i?(p, q, k) - fc'^AJp)i?(p, q, fc) + • • • 

- i?(p, q, fc) + fc'^AJp)i?(p, q, fc) + fc^^9bA„(p)i?(p, q, fc) 
= - k'^d.MP, Q, k) + k-q'dMp)Rip, q, fc) + • • • 



Rip, q, fc) 1 - k'^A^ip) - k'^q'dMp 



k" 

q 



'k' (daT^f) ip) + -q'q'k^'' {daTSL) iP) + ■■■]+ k^'dMp) [q^k= ) (p) 

(d.T^f) (p) + q^q^k"^k- {daT§^ (p) + q'kWk^dMp) (x^f) (p) + • • • 



(E9) 
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where the summation convention over the repeated indices a, 6, z, j, to = 1, • • • ,d is imphed. 



At last, we are in a position to evaluate the terms contributing to the function in Eq. (El I. We start with 



fl.k ^IYI ^q(i^P-k>Sp^q¥Pq,SPk] 



P' 



q Xp' Xp'+/c+g 



(ElO) 



q p.p' 



q p,p' 



i E E [ ^ 9) • ^(P)] [ A q) • B{p' 



q p,p' 



Xp Xp— fc — q Xp/ Xp'-|-fc+q 



Xp \p—k — q Xp' Xp' + fc+q 



where we used that = e^^"^djA^ = \e^^™-Fj^ or, equivalently, F^^ = e^^^i?™. We now break the Berry field 
strength into two contributions, i.e., B{p) — B + SB{p). If so. 



fi,k = 5] 5] t^q [ (fc A q) • (B + SBip)) ] [ (fe A q) • (B + 

"3 P,P' 

= E E ^^[(fc^?)-^]!^^^?)-^ 

q p,p' 

EE«<3[('^^9)-^][(^^9)-^^(p') 

q p.p' 

EE«<3[('^^9)-^^(p)][('^^9)-^^ 

q p.p' 

- ^ E E [ ^ 9) ■ ^^(p)] [ (fe ^ 9) ■ ^B{p') 

q p.p' 

= -^E"«[(^^'?)-^][(^^'?)-^ 

-^EE«'j[(^^'?)-S][(fcAg).5B(p'; 

q p' 



Xp Xp— fc— g Xp' Xp'+fc+q 



Xp y^p—k — q Xp' Xp'+fc+q 



Xp Xp— fc— q Xp' Xp'+fc+q 



Xp )Cp—k — q Xp Xp'+fc+g 



Xp \p—k — q yip' Xp' + fc+q 



(Ell) 



P-k-q Pk^q 



P—k — q Xp' Xp'+fc+qf 



Xp Xp— fc — g Pk-\-q 



q p,p' 



q p^p' 



Xp \p—k — q Xp' Xp'+fc+q 



In a uniform liquid-like ground state we have (p^ ) oc 5^ q and, due to the relation fc'^g'' i^^^^ = {k A q) ■ B, we can 



replace p±k±q by As a consequence, we can drop the first three terms on the last equality of (Ell I up to 

order q^k^. We are then left with: 



fi,k = E E ^.J [ ^ 9) • SBip)] [ (fc A q) • SBip' 



q p,p' 



(E12) 
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where fip = Xp Xp is the number operator projected on the lowest band. Similarly, 

^\^^<1 {\^Pk:5p-q][5pq:5P-, 



while 



— /i,-fe 

= -^EE«'j[(^^'?)-'5i?(p) 



{kAq)-SB{p') 



q p,p' 



fa.k = ^ E ( ^P-Q ^^P-k' t^^q' ^Pk] ] 



p+q 



(E13) 



(E14) 



The matrix element 6p_qXpXp+q^ vanishes in the limit q -> and, therefore, the only term that contributes to 
/s k up to order q^k^ is 



— q Xp Xp+q 



q p 



q p 



1 r / \ 1 / 

2 E E ^q [(^^9) • (a^J iP)\k''{Sp-qXlXp+q 



The condition (E5l implies that f^j^ — f^j^^ which then dehvers 

1 r / OB \ 1 / 

/4,fc - 2 E E ^'q [ - ^ 9) • ( ) J { ^l+q %^Pq 



q p 

where we have used that {B{p))* — —B{p). 

Putting together all the contributions, we obtain 

/a. = - E E "q [ ^ 9) • SB{p)] [ (fe A q) • SBip' 
q p,p' 



EE«q 



q p 



dB 



if' ^ Q) ■ [ dp^ ) ^p ^p+q) " ^ • ( TwT ) \^p+q ^p "^^q 



dB 



dp" 



(E15) 



(E16) 



(E17) 



where the summation convention over the repeated indices a = 1, • • • , is implied. Finally, the analytical continuation 
B = \B delivers Eq. ( |4.10a[ ). 
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